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We define a family of binary outcome n-party m < n settings per party Bell inequalities whose members 
require the least detection efficiency for their violation among all known inequalities of the same type. This 
gives upper bounds for the minimum value of the critical efficiency — below which no violation is possible — 
achievable for such inequalities. For m = 2, our family reduces to the one given by Larsson and Semitecolos 
in 2001. For m > 2, a gap remains between these bounds and the best lower bounds. The violating state near 
the threshold efficiency always approaches a product state of n qubits. 
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I. INTRODUCTION 

In a Bell experiment distant parties perform measure¬ 
ments on a shared physical system. From very reasonable 
assumptions it follows that the correlations between their 
results must satisfy certain inequalities. Quantum mechan¬ 
ics predicts that often these Bell inequalities will not be re¬ 
spected |[l],|2] ■ Their violation, which excludes all local real¬ 
istic explanations for the quantum world, and which is often 
referred to as nonlocality, is one of the most surprising and 
counterintuitive features of quantum physics. One of the 
earliest experimental demonstrations of the Bell violation 
were done by Aspect et al. lilt], and many experiments on 
various systems have been performed since (see, e.g. Ref. 
101). Unfortunately, due to technical imperfections there 
are still loopholes, which in principle allow local realistic 
explanation for the experimental results. 

An important loophole, the so-called detection loop¬ 
hole Mpears in experiments performed on systems of pho¬ 
tons 1^. Some photons get lost during transmission and the 
efficiency of the present day photon detectors is also lim¬ 
ited. If the proportion of photons detected falls below a cer¬ 
tain critical value r]crit, which depends both on the Bell in¬ 
equality and on the system, then the results become compat¬ 
ible with local realism. As an important property character¬ 
izing the inequality and the system, Vallone et al. 11 intro¬ 
duced the notion of the robustness of nonlocality R. They 
defined this quantity as the maximum fraction of detection 
events that can be lost such that the remaining ones still do 
not admit a local model. It is easy to see that R= 1 — rjcrit- 
In the present paper it will be assumed that the detection ef¬ 
ficiency is the same for all parties. There are experiments, 
for example those involving entangled photon-ion systems, 
where some of the parties (the ones working with the ions) 
have virtually no loss of events. Then the experiment tol¬ 
erates higher losses for the remaining parties than in the 
symmetric case considered here H. 

Another important loophole is the locality loophole, 
whose closure requires space-like separation between the 


parties. This can realistically be ensured only in Bell tests 
performed on photon systems. Each loophole has been 
closed by some experiment, but not all simultaneously, for 
example, no single experiment closed the locality and the 
detection loophole at the same time (for a comprehensive 
review, see Ref. |0]). Such a verification would be interest¬ 
ing and fairly important from a philosophical point of view, 
although — as so many predictions of quantum mechanics 
has proven correct — at present virtually everyone is quite 
convinced that the Bell violations observed so far are gen¬ 
uine, and not results of some conspiracy plotted by nature, 
using different loopholes in different situations to mislead 
us. 

However, there are also practical reasons why loophole 
free Bell violation would be important to achieve. Bell vi¬ 
olation is the basis of the so-called device independent pro¬ 
tocols of mantum information technology (see for exam¬ 
ple Ref. 121 for a recent review). Such protocols, based 
on the fact that Bell violation is very hard to fake, would 
make it possible for anyone to check the proper functioning 
of an apparatus without knowing much about the details of 
its internal structure. Therefore, construction faults would 
be easier to notice and it would be harder for a malicious 
manufacturer to build back doors into the device. This ap¬ 
proach promises to perform cryptographic tasks with an un¬ 
precedented security d, produce genuinely random num¬ 
bers fH , and carry out black-box state tomography ifl^ . 

The violation of the simplest Bell inequality, the two- 
party two-setting per party (n = m = 2) Clauser-Horne- 
Shimony-Holt (CHSH) inequality 1101 requires at least 
Rcrit = 2/3 ~ 0.667 detection efficiency 13. This 
value may be reached with a partially entangled two qubit 
state, which approaches a product state near the thresh¬ 
old efficiency. In this case the robustness of nonlocality is 
R = 1/3. As far as we know, with two parties perform¬ 
ing measurements on a pair of qubits no better value has 
been achieved. In Ref. 13 ^ bipartite four-setting inequal¬ 
ity was found with a somewhat lower critical efficiency of 
rjcrit = ( 3/5 — l)/2 ft! 0.618. To get this value a pair of 
ququarts had to be used. Larsson and Semitecolos Il3 have 
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introduced a family of n > 2 party binary Bell inequalities 
with m = 2 settings per party with pcrtt = n({2n — 1). 
They have also proved that this value is optimal, no two- 
setting n-party Bell inequality may be violated with any 
lower detection efficiency. This means that for m = 2 an 
efficiency of less than 1/2 is never enough. 

Let p* be the smallest of the critical detection efficiencies 
belonging to the Bell inequalities characterized by a certain 
number of parties and numbers of settings and measurement 
outcomes for each party. Above ry* there are quantum mea¬ 
surements exhibiting correlations that can not be explained 
by any local hidden variables model. The value for rjcrit 
given by Ref. d] and cited above is just p* for binary- 
outcome n-party two-setting per party Bell inequalities. For 
n-party and more than two-setting per party (m > 2) in¬ 
equalities 77* is not known. Massar and Pironio d has 
given a lower bound for this quantity by constructing ex¬ 
plicit local hidden variables models reproducing the cor¬ 
relations for efficiencies below this bound. For m = 2 
the bound is the same as the one of Ref. fl^. B ased on 
combinatorial considerations Buhrman et al. fl^ has de¬ 
rived an upper bound for 77*, when m is some power of 
two. Ref. lU^ also provided upper bounds for many (n,m) 
combinations (m prime) by giving explicit Bell inequalities 
corresponding to those efficiency values. The inequalities 
of that work may become useful in practice as well, because 
they require the Greenberger-Horne-Zeilinger (GHZ) state 
1 ^, very well distinguishable measurement settings even 
near the threshold efficiency, and the noise tolerance is also 
quite good. 

Very recently tripartite Bell inequalities have been con¬ 
structed requiring low detection efficiencies for their 
violation with the W state ll23l . With three measurement 
settings for one of the parties and two settings for the other 
two parties a critical efficiency of rjcrit = 0.6 has been 
achieved. The same efficiency has been necessary for an 
inequality symmetric for the permutations of the parties 
with 777 = 3 settings per party. For m = 4, 6 and 8 this 
value has been improved to rjcrit = 0.509036, 0.502417 
and 0.501338, respectively. For larger m the construction 
would have required too much computational resources but 
the result would probably have remained above 0.5. How¬ 
ever, if a small admixture of a product state to the W state 
has been allowed, rjcrit = 0.5 has been found for an m = 4 
inequality. This value could not be improved by increasing 
m. 

In this work, we have used an iterative procedure alter¬ 
nating a linear programming step and a semidefinite pro¬ 
gramming one to get the tti = 3, n = 3 and n = 4 Bell 
inequalities with the smallest possible rjcrit- In this pro¬ 
cedure we made assumption neither about the symmetry 
of the inequality nor about the properties of the state, not 
even about the dimensionality of the Hilbert space. We have 
got tripartite inequalities reaching rjcrit = 0.5 already with 


three settings per party, and for tti = 3, n = 4 we have ob¬ 
tained rjcrit = 6/13 ~ 0.46154. We could get the same val¬ 
ues with confining ourselves to permutationally symmetric 
inequalities. Numerical calculations showed that the state 
giving the maximum violation near the threshold efficiency 
is a state of n qubits, namely a mixture of the W state and 
a product state that approaches the product state near rjcrit, 
the same behavior as the one observed for the two-setting 
inequalities of Larsson and Semitecolos in. The measure¬ 
ment settings could be chosen the same for all parties. Each 
measurement operator could be characterized by a single 
angle that approached zero at rjcrit, but at different paces 
for the different settings, again similarly to what is seen in 
Ref. iS. These observations helped us to make the gen¬ 
eralization and to define a whole family of Bell inequalities 
with any m and n > m. The members of the family require 
the least detection efficiency for their violation among all 
known inequalities of the same type, giving upper bounds 
for 77*. For m = 2 the inequalities are the same (apart from 
swap ping of the two measurement settings) as the ones in 
Ref. in. 

The structure of this paper is as follows. In Section HU 
we give explicitly the family of multipartite Bell expres¬ 
sions. Then, in Section uni we give the optimal quantum 
violation of these Bell inequalities with qubit systems ac¬ 
counting for finite detection efficiencies. Our main result is 
presented in formula ® providing (to the best of our knowl¬ 
edge) the best upper bound on 77* for any n > m parties 
and m settings per party. To get this formula, we have used 
two claims; one about the equality of certain quantum con¬ 
ditional probabilities and the negligibility of others, while 
the other concerning the classical bound. The first claim is 
proven in Section |I3 while the second one in Section |V] 
To demonstrate the behavior of the relevant quantities, we 
highlight the example 77 = tti = 3 in section [VT] Here, we 
give the explicit Bell expression and also show the optimal 
measurement settings (angles) and Bell violation departing 
from the value of the critical detection efficiency. Finally, 
we summarize our conclusions. 


II. THE INEQUALITIES 

We consider the Bell scenario with n distant observers. 
Each observer may freely choose between a set of m binary 
measurement settings with possible outcomes zero and one. 
Let Aij denote the jth measurement of the zth party. A Bell 
inequality corresponding to this setup may be written as; 

mm m 

EE-E ■ ■ ■ r^njn) — 

jl=0j2=0 jn=0 

( 1 ) 

where L is the classical bound, and 
■ ■ ■ -^njn) = 7^(11 . . . I|^ljl^2j2 ■ ■ ■ -^njrt) 
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denotes the conditional probability of all parties getting 
outcome one given they have measured Aij^, ^ 2 ^ 2 ’ ^^id 
so on. For the sake of brevity of the notation, for terms 
involving only a subset of observers, we have introduced 
a zeroth measurement whose outcome is always one, 
and which is always supposed to be performed. This 
measurement appears in these terms for the observers not 
members of the subset considered. 

The inequalities we introduce in the present paper are in¬ 
variant with respect to the permutations of the parties, and 
n > m holds. Besides terms involving all n observers, 
they contain only terms with n — 1 parties (that is terms 
containing equal or less than n — 2 parties are zero). Due 
to the permutational symmetry. Bell coefficients having the 
same set of indices in different orders have the same value. 
Therefore, we only give explicit values for coefficients with 
indices ji < j 2 < js • • • < in- 

The nonzero (n— l)-party Bell coefficients have the value 
of minus one. Besides the index zero for the party not in¬ 
volved, they have indices from 2 to m once, plus index m 
another n — m times if n > m. For example, if m = 2 
and n > 2 then B 02 . . .2 = —1, or if m = 3 and n > 3 
then i? 023...3 = —1- For larger values of m the nonzero 
coefficients are 

F^0234...mm...m — 1; (2) 

and the ones we get by permuting the indices of the coef¬ 
ficients above. Thus there are n!/(n — m -I- 1)! nonzero 
(n — l)-party coefficients. 

Furthermore, we have n-party Bell coefficients whose 
value is positive. They are the ones that have the same set of 
indices as the nonzero (n — 1) -party coefficients above, but 
index zero is replaced by anything from one to m. If this 
index is m, the value of the coefficient is (n — m -I- 1), oth¬ 
erwise it is one. Therefore, we may write (again for larger 
m) 

— 1 

B2234.. — 1 

B2334.. — 1 

B2344.. — 1 

f^234...mmm...m = ri 777.-1- 1. (3) 

Again, all coefficients we get by permuting the indices have 
the same value. Thus the first and the last line represent 
77!/(77 — 777 -I- 1)! and 771/(77 — 777 -I- 2)! Coefficients, respec¬ 
tively, while the lines in betweenrepresent 771/2(77 — 777 -|-1)! 
coefficients. 

There are also some 77-party Bell coefficients having neg¬ 
ative values; 

(4) 


where the indices are such that the quantum values of the 
conditional probabilities multiplying these coefficients cal¬ 
culated with the quantum state and the measurement oper¬ 
ators we will give later are vanishingly small compared to 
the quantum conditional probabilities multiplying the coef¬ 
ficients given by Eqs. da and da. Then the actual values 
of these coefficients will not influence the quantum value of 
the Bell expression (the left hand side of Eq. ([TJ, therefore 
we may choose them freely. As we will show later, this free¬ 
dom makes it possible to achieve that the classical bound L 
appearing in Eq. dD is zero, which is appropriate for our 
purpose. There are many solutions, and they all give the 
same value for the threshold detector efficiency. The maxi¬ 
mum quantum violation above the threshold efficiency does 
depend on the actual choice, but we have made no attempt 
to find the best one. 

All other Bell coefficients are zero. 


III. QUANTUM VIOLATION WITH DETECTORS OF 
LIMITED EFFICIENCY 

Let us consider the quantum violation of our Bell in¬ 
equality. If we confine ourselves to von Neumann measure¬ 
ments on a pure ?7-qubit state It/), the quantum value of the 
conditional probability appearing in Eq. dD may be written 
as: 

n 

F’Q (^ 171 ^ 2^2 • ■ ■ -^njn) ~ (V^l 

i=l 

where Aij for 1 < j < 777 is the measurement operator 
corresponding to measurement Aij, = / is the identity 
operator in the subspace of the observer, and Ay . de¬ 
notes the tensor product of the operators of the parties. The 
state 1 ^) we consider in this paper is invariant with respect 
to the permutations of the parties, and all parties have the 
same set of measurement operators, that is Ay = Aj is in¬ 
dependent of i. Then, if the Bell inequality of Eq. dD is 
permutationally invariant and the classical bound is zero, 
the condition for its quantum violation can be written as: 

mm m 

E E-- E J2; ■ • ■ I jTl)f?7i72. .7n 

71 = 072=71 jn=jn-l 

n 

('01 AjAV') > 0 , (6) 

i=l 

where 7r(ji , j 2 , ■ ■ ■, jn) denotes the number of independent 
permutations of its arguments. Let of suppose that the ob¬ 
servers have detectors of limited efficiency 77 < 1. Let they 
agree that each of them signal outcome zero whenever the 
particle is not detected. In this case each term in the left 
hand side of Eq. dH) must be multiplied by the probability 


Bk ik2k3k4,...kn ^ 0 ; 
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of detecting all particles concerned to get the condition for 
the detection of the violation. This probability is p* for an l- 
party term, where / is the number of nonzero indices. At the 
threshold efficiency rjcrit the resulting expression is zero. 
This condition gives an nth order equation for pcrit- In our 
case the solution is simple, because we have restricted our¬ 
selves to n-party and (n— l)-party terms. The condition for 
being able to detect the violation is 

mm m 

32 = 1 33 =32 3n =3n-l 
n 

i=2 

mm m 

E-.E J2l ■ ■ ■ ! •■■jn 

il=lj2=il jn=3n.-l 

n 

(v^i(8)4i^) >0. (7) 

i=l 


IV. THE QUANTUM CONDITIONAL PROBABILITIES 


Let all measurements be real ones, that is performed in 
the X — Z plane. Then each operator Aj can be character¬ 
ized by a single real variable (pj: 


^■|0) 


^.|1) 




sin tpj 

1 -I- cos ^ 
2 


|l) = c-|0) + s,|l) 
^|l) = s,|0)+c+|l). 


(9) 


A zero angle corresponds to the measurement giving out¬ 
come one with probability one for the 11) state. Let all mea¬ 
surement angles be small, and let they obey the following 
hierarchy: 

0 < < I02| < |«!>3| < ••• < (10) 

Let the state be 


Then rjcrit is minus one times the ratio of the (n — 1)-party 
and n-party contributions in the expression above. The for¬ 
mer must be negative while the latter positive such that the 
left hand side is positive above rjcrit- 

In the next section we will show that with an appropriate 
choice of the state vector and the measurement settings, the 
quantum conditional probabilities associated with each Bell 
coefficient given in Eqs. (|2l) and Q are equal. As the same 
number appears both in the numerator and the denomina¬ 
tor of the ratio giving ricrit, we can simplify with it. The 
factor TT giving the number of independent permutations of 
the indices is n!/(n — to -b 1)! for the {n — l)-party coef¬ 
ficient and for the n-party coefficient given in the first line 
ofEq. dill. The factor is one half of this for the next to — 2 
coefficients, and it is n!/(n — to -b 2)! for the last one. Tak¬ 
ing into account the values of the coefficients we can give 
the expression for rjcrit as: 


^crit — 


1 _1_ m—2 
^ 2 


n—m+1 
n—m+2 


2 -b TO — 


( 8 ) 


n—m+2 


|'0) = cosq;|0) — sina|lU), (11) 

where 

| 0 ) = | 00 ... 0 ) ( 12 ) 
\W) = E(|10...0) -b I01...0) -b-- - -b I00 ...1)). 

Vn 

(13) 


This state is permutationally symmetric. The angle a will 
be small, therefore the state is predominantly the separable 
|0) state with a very small amount of \ W) state added. 

Eirst we calculate the quantum conditional probability 
multiplying the (n — 1)-party Bell coefficient given by 
Eq. (E). Erom Eqs. (E and (fTTli we get: 

Pq{AioA 222^33 . . . ■ ■ ■ A^m) — 

cos^ aSoo — 2 COS a sin aSwo + sin^ aSww, (14) 


Here we have simplified the fraction by the common fac¬ 
tor n!/(n — TO -b 1)!. This expression is our main result. 
Our Bell inequalities are such that their violation, at least 
in principle, can be detected by using detectors of as low 
efficiency as given above. Therefore, this is an upper bound 
for the threshold detector efficiency which may be achieved 
for n-party m <n setting binary outcome Bell inequalities. 
The expression gives 2 /(to -b 1) for n = m and 2 /(to -b 2) 
in the limit of large n. However, if we increase n, we can 
achieve less than by increasing both to and n by just one. 

To prove Eq. (l8]l we must show that the relevant quan¬ 
tum conditional probabilities are really equal as we have 
claimed, and that the local bound is zero. We will do this in 
the next two sections. 


where 

n 

= (ctIJ® (^ + |t), (15) 

i=2 

with Ai = Am when i > to. As the first step to calculate 
these matrix elements, from Eqs. (E and (fT2l i we get: 

n n 

/®(g) + |0) = |0)®(g)(c-|0)+s+)). (16) 

i—2 i—2 

Then, by using Eqs. (fTSI) and dnii it is easy to calculate 
^oo and S'vvo- We may also get Sww from a similar, but 
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somewhat lengthier calculation. The result is: 

n 

Soo = n 

i=2 

^ n n ^ 

-L TT - *^00 

"T 


Sww = So 


k=2 fc i=2 

1 1 

- + 7^ 
n 


(17) 


As |cj I <C I Si |, ^ is a small number. By substituting 
Eqs. (fTTl i into Eq. (fT4li we get 


(^10^22^33 ■ • ■ ■ • ■ -^nm) — 


5, 


00 


^ 2 -2 
Sin a \ sin a 


cos a — 


e 


(18) 


The first term in the square bracket dominates over the sec¬ 
ond one. However, as we will see later, the contribution 
from the n-party terms is of the same order as the second 
term. The threshold efficiency may only be finite if the first 
term vanishes. We can achieve that by choosing the mixing 
angle a characterizing the quantum state such that tga = 
With this choice 

TD ! A A A A A A \ ® 

(.-^ 10 -^ 22-^33 • ■ • ■ • ■ -^nm) — • 

(19) 

Now let us calculate a general n-party conditional prob¬ 
ability. Analogously to Eqs. (fT4l i and (fTST i we may write: 

-P(3(^ljl^2j2 ■ • ■ ^njn) = 

cos^ aFoo — 2 cos a sin aFwo + sin^ aFww, (20) 

where 


Far — (o’I Aij. |r). 

i=l 


( 21 ) 


We can calculate these matrix elements similarly to the ones 
given in Eq. (Ell, and we get; 


"00 


n< 


n n 

Sfc i=i ^ 


Fww = 




( 22 ) 


By substituting Eqs. (l22li into Eq. (l20l i we get 


7^q(Aij^A2j2 ■ • ■ — Fqo I COSO; 


sin a 


Fqo sin^ o [ ^ - 1 • (23) 


C X 


Here we used the relation tgo = 

Let us consider the quantum conditional probabilities 
multiplying the n-party coefficients whose values are given 
explicitly in Eq. ([3). The product giving Fqo includes ex¬ 
actly the same elements as the one giving ^oo, plus one 
more (see Eqs. (fTTl i and (l22l i. The same is true for the 
sums defining x und The extra element is the one of 
index I for the Zth line of Eq. (13- Then Fqo = SqqcJ' and 
1/x = 1/f + si/{cY a/u). By substituting these values into 
Eq. (l23l l we get that the result is jcY times the value given 
in Eq. (fTSl l for the (n— 1)-party conditional probability. But 
this extra factor is one in the limit of small angles. There¬ 
fore, for our state and measurement settings, the quantum 
conditional probabilities multiplying each of the Bell coef¬ 
ficients of Eqs. (13 and ([3 are equal indeed, and have the 
value given in Eq. ( fT9] l. 

Now, let us consider which are the Bell coefficients 
(Eq.|4|i whose value may freely be chosen. Erom Eqs. ( fT9b 
and (123 it follows that the quantum probability character¬ 
ized by indices ki < k 2 <■■■< kn is negligible compared 
to the significant ones if 


n n — / n 

(V — 

nr= 2 c- \Uc- 



< 1. 


(24) 


Here we substituted the values of Soo, Fqo, C und x from 
Eqs. (Ell and (123 . However, c“ = (1 — cos(^i)/2 r; 
and Si = —sin0^/2 r; if the angles are small, 

therefore the above condition may be written as: 


nr=2 



(25) 


If ki = 1, due to the hierarchy imposed by Eq. (fTOl i the 
dominant term in the parentheses is \/4>i. Then the condi¬ 
tion above may be written as (nr =2 «1- 

This relation holds if tpki < (j^i for 2 < i < n, and for 
at least one index w the inequality is strict, in which case 
4>k^ ^ (l>w- When deriving the formulae we have implic¬ 
itly supposed that (j)i ^ 0. However, the result is true for 
= 0 as well. lfki = i + \ for i < I, while ki = fc;_i = /, 
the dominant term in the parentheses is 1/0; and the condi¬ 
tion becomes (nr=i-i-i <g; 1. Then if the 

same relations hold between 0i and 0^^ as for the ki = I 
case, but now for I + 1 < i < n, Eq. (123 will hold and the 
threshold efficiency will be independent of the choice of the 
associated Bell coefficient. 
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V. THE CLASSICAL CASE 

For an n-party Bell inequalitiy with m measurement set¬ 
tings per party the classical bound is zero if 

mm m 

EE-E - • ■ • ^njn ^ 0) (26) 

Jl =0^2=0 jn=0 

for all deterministic strategies, which are defined by the ac¬ 
tual choices of (i = 1,..., n; j = 1,..., m), where 
each of them may take the value of zero or one, while 
flio = 1 for all i. The a^o appear in the subcorrelation 
terms. A deterministic strategy means that each party has 
a definite outcome for each measurement setting with prob¬ 
ability one. This outcome for the jth measurement of the ith 
party is given by . Let us denote the matrix with elements 
Uij by a. Then for permutationally symmetric inequalities 
Eq. (l26l l may be rewritten as: 


a factor of 2(n —m-f 1)!. The (5(d) > 0 denotes the contri¬ 
bution of terms due to the negative n-party Bell coefficients 
of Eq. (HI, multiplied by minus 2{n — m + 1)1. 

A permanent is independent of the order of the rows 
and columns of the matrix. Therefore, may rearrange the 
columns of the matrices appearing in Eq. ( l29b such that 
their jth column is the (j -I- l)th column of d if j < m — 1, 
the rnth column ifm — 1 <j < n — 1 and only their 
last column differ from each other. However, the sum of 
two permanents whose matrices differ from each other only 
in one single column (row) is the permanent of the matrix 
having the sum of the those columns (rows) in that position. 
Eurthermore, multiplying a permanent with a factor has the 
same effect as multiplying one of the columns (rows) of the 
matrix. Therefore, Eq. ( |29] | may be rewritten as: 

perm(^(d)) — Q{a) < 0, (30) 

where 


32 33 


EE-E Ji, J2) ■ • ■ jjn) ^ 0, 

31=0 32=0 3n.=0 


where 

^(jl I J2 7 • ■ • : Jn) 


^(jl 7 J 27 ■ • ■ 7 Jn) 


n! 


O-eSn i=l 

perm(a(d, ji,j27 ■ ■ ■ 7 jn))- 


(27) 

n 

E ^ 

(28) 


The sum in the first line extends over all elements a of the 
symmetric group S'„, that is over all permutations of the 
numbers 1 ,2,...n. We should only have summed over the 
independent permutations of the numbers ji, j 2 ■ • ■ jn- The 
difference is taken care of by the factor in front of the sum¬ 
mation with 7r(ji, ^2 7 • • ■ 7 jn) denoting the number of in¬ 
dependent permutations of its arguments. The sum itself 
is nothing else but the permanent of the n x n square ma¬ 
trix a{a,ji,j 2 , ■ ■ ■ ,jn) consisting of the jith, j 2 th,..., j„th 
columns of the d strategy matrix. 

Eor the Bell inequality given by Eqs. (|2]i-(|4li the condition 
Eq. (l27l) that the classical bound is zero may be written as: 


2perm(Q;(d, 1234 ... mm ... m))+ 
perm(Q;(d, 2234 ... mm ... m)) + 
perm(a(d, 2334 ... mm ... m)) + 
perm(Q;(d, 2344 ... mm ... m)) -f • • • -f 

_ jYl —j— 1_ 

2 --perm(Q;(d, 234 ... mmm ... m))— 

Th TTZ H” ^ 

2perm(Q;(d, 0234 ... mm ... m)) — Q{a) < 0. (29) 

Here we have used Eq. (l28l l. we have substituted the actual 
values of Tr{ji,j 2 , ■ ■ ■ ,jn) and simplified the inequality by 


Pij 

Pij 


— ^im 


j < m — 1 
m—l<j<n—1 


m —1 

= 2aii + Qik 

k=2 


n — m + 1 

I --( 

n — m + 2 


(31) 


Our Bell inequality is appropriate if Eq. (I29]) holds for 
all deterministic strategies. The permanents appearing in 
the expression can not be negative, as all the elements of 
the matrices involved are zero or one. Therefore, strategies 
giving zero for all terms with a positive sign trivially sat¬ 
isfy Eq. ( |29] |. We also do not have to consider explicitly 
any strategy leading to a nonzero factor for a Bell coeffi¬ 
cient whose value does not affect the threshold efficiency 
(Eq.|4]i. The Q{a) depends linearly on such a coefficient, if 
we choose its absolute value large enough, we can always 
satisfy Eq. (l29l l. 

Eirst let us consider strategies giving a positive value for 
the first term of Eq. (l29l l. This is true if and only if there is 
at least one party having outcome one for the first measure¬ 
ment setting, at least one other party having outcome one 
for the second measurement setting, and so on, and finally, 
there are at least n—m+1 other parties having outcome one 
for their mth measurement setting. Eor a permutationally 
invariant inequality all strategies differing from each other 
only in the order of the parties give the same constraints for 
the Bell coefficients. Therefore, without sacrificing gen¬ 
erality, we may assume that the first party’s first measure¬ 
ment outcome is one, the second party’s second measure¬ 
ment outcome is one, and so on, and finally, the last n — m 
party’s mth measurement outcome is one, that is Uik = 1 
where k = min(i, m). We only have to consider strategies 
whose matrices contain zeros to the left of these elements 
in each row, that is aik' = 0 if 0 < k' < k. If Uik' were 
one, then the factor of Bi 2 ,,,(i-i)k 'would not 
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be zero. However, the quantum conditional probability as¬ 
sociated with this coefficient is {4>k'/4>k)‘^ ^ 1 times the 
one associated with the significant coefficients (this follows 
from Eq. (l25T l. therefore, we need not consider such a strat¬ 
egy explicitly indeed. The matrix of the remaining strate¬ 
gies for TO = 6 and n = 9, without the zeroth column is: 


/I 

ai2 

ai3 

ai4 

Ol5 

016 \ 

0 

1 

023 

024 

O25 

026 

0 

0 

1 

034 

O35 

036 

0 

0 

0 

1 

O45 

046 

0 

0 

0 

0 

1 

056 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

Vo 

0 

0 

0 

0 

1 / 


With such a strategy matrix (5(a) appearing in Eqs. ( |29] | and 
(l30t is zero: terms contributing to it have zero factors. It is 
also clear that increasing the number of parties beyond to 
does not allow any more free parameters. 

Eor this example Eq. (l30l) may be written as: 



/ 

012 

Ol3 

014 

015 

016 

016 

016 

016 

Pin \ 



1 

O23 

024 

025 

026 

026 

026 

026 

P2n 



0 

1 

034 

035 

036 

036 

036 

036 

P^n 



0 

0 

1 

045 

046 

046 

046 

046 

Pau 

perm 


0 

0 

0 

1 

056 

056 

056 

056 

Pbn 



0 

0 

0 

0 

1 

1 

1 

1 

Pmn 



0 

0 

0 

0 

1 

1 

1 

1 

Pmn 



0 

0 

0 

0 

1 

1 

1 

1 

Pmn 


V 

0 

0 

0 

0 

1 

1 

1 

1 

Pmn j 


(33) 

Let us expand the permanent along its last row. There are 
n — m + I ones in that row, all in columns equal to each 
other. These terms of the expansion give n — to -I- 1 times 
the permanent of the matrix with that colunm and the last 
row deleted. The contribution from the last element will 
be Pmn times the permanent of the matrix with its last row 
and last column deleted. The sum of the contributions may 
be expressed as the permanent of a single matrix whose ith 
element in its last column is {n — m + l)^i„ -|- /3mn<iim, 
which we can calculate from Eqs. (OTI) and (l32]) . The value 
we get is 


C m—1 
2aii + ^ ^ (lik 2 
k^2 

(34) 

which is nothing else that n — m -I- 1 times the last column 
of the matrix with n — 1 parties. Therefore, the value of the 
permanent for n parties is n—m +1 times the value for n—1 
parties. If we repeat this n — m times we get (n — to -I- 1)! 
times the permanent for n = m. 

Let us consider n = m. Eor to = 6 Eq. (l30l) may be 


n — m + 1 


- 2 


written as: 



/ 

012 

Ol3 

Ol4 

Ol5 

O16 

E.2 

Oli 


\ 




1 

O23 

024 

O25 

026 


02i — 

1 



perm 


0 

1 

O34 

O35 

036 

nU 

03i — 

1 


< 0 



0 

0 

1 

O45 

046 

045 -f 

046 — 

1 




1 

0 

0 

0 

1 

056 

056 — 

1 


) 



0 

0 

0 

0 

1 

-1 





(35) 

The values in the last column follow from Eq. (OTI) . If we 
expand the permanent along its last row we get the perma¬ 
nent of a matrix of the same form corresponding to to — 1 
measurement settings. If we repeat this m — 1 times, fi¬ 
nally we get zero for the permanent. Therefore, Eq. (OOl) and 
equivalently, Eq. (l29l) is satisfied as an equality. Although 
we demonstrated this result on the example of m = 6 , it is 
obvious that it is true for any to. 

So far we have only dealt with strategies that lead to a 
nonzero factor for the first term of Eq. ( l29b . Now let us 
consider strategies giving zero for the first I — 1 terms of 
Eq. (l29l l and nonzero for the (th term. If n = to = 8 and I = 
5 the strategy matrix to be considered without the zeroth 
column is the following: 


( ail 

1 

Ol3 

Ol4 

Ol5 

016 

Ol7 

O18 

021 

022 

1 

024 

O25 

026 

O27 

028 

031 

032 

O33 

1 

O35 

036 

O37 

038 

0 

0 

0 

0 

1 

046 

O47 

048 

0 

0 

0 

0 

1 

056 

O57 

058 

0 

0 

0 

0 

0 

1 

067 

068 

0 

0 

0 

0 

0 

0 

1 

078 

V 0 

0 

0 

0 

0 

0 

0 

1 


Elements = 0 and aki = 0 for k < I, otherwise 

the fcth term of Eq. (l29l l would be positive. Eurthermore, 
Qij = 0 if j < i for i > /. Eor example, 073 above is zero. 
If it were not, the factor multiplying 7723345568 would be 
positive. The quantum conditional probability associated to 
this coefficient is small by a factor of <C 1 (see 

Eq. (|25 T i). Eor the interesting strategies a larger number of 
parties would mean an additional n — m copies of the last 
line, which contains a single one in the last position. Only 
this property of the matrix was used when we showed that 
increasing n beyond to only introduces a positive factor, 
which does not affect the validity of the inequality. There¬ 
fore, it is enough to deal with n = m. 











For the present example Eq. (EOll may be written as: 


/ 1 013 ai4 015 Oi6 017 0,18 Plm \ 

022 1 024 025 026 O27 O28 /32m 

O32 O33 1 O35 O36 O37 O38 /33m 

0 0 0 1 046 O 47 048 Pim , n 

perm n o o i a ^0. 

0 0 U 1 056 O57 058 P 5 m 

0 0 0 0 1 067 068 /36 m I 

000001 078 /37m 
V0000001 /38m/ 

(37) 

This permanent is equal to the product of the permanents of 
two blocks along the diagonal of sizes (/ — 2) x {I — 2) and 
{n — I + 2) X (n — / + 2). In the present case 


perm 






/I 

046 

047 

CO 

^4m \ 

(1 

Ol3 

014 

\ 

1 

056 

057 

058 


022 

1 

024 

xperm 

0 

1 

067 

068 

i^Am 

\ 032 

O33 

1 7 

/ 

0 

0 

1 

078 

i^Am 





\0 

0 

0 

1 

^Am / 


(38) 


If we substitute the values of the last column of the second 
matrix from Eq. OTl i. we get a matrix whose structure is 
very similar to the one in Eq. dTSl) . There is one important 
difference. Now the last elements of each row is always the 
sum of the preceding elements minus two. In the case of 
Eq. drsl) this is not true for the first row, the minus two is 
missing. Due to this difference, if we follow the steps we 
have taken in the case of Eq. (EB, we will get not zero, but 
minus two. Therefore, Eq. ( |29] | is satisfied now as a strict 
inequality. 


VI. THE 333 CASE 


Eor m = n = 3 a permutationally invariant inequality 
with the properties required may be defined by the follow¬ 
ing coefficients: 

Bo23 = —1 Bii2 = —1 Bii3 = —1 Bi22 = —2 

Bi23 = +1 B223 = +1 B233 = +1. (39) 

Bell coefficients having the same set of indices in different 
orders as any of the coefficients above have the same value. 
All other coefficients are zero. 

We have got this inequality by iterating a linear program¬ 
ming step and a semidefinite programming one. The for¬ 
mer one provides the Bell inequality of zero classical bound 
with no single-party marginals having the lowest critical de¬ 
tector efficiency from known matrix elements of the mea¬ 
surement operators. The procedure is the same as the one 
we used in Ref. . In the semidefinite programming step 
we applied the method of Navascues, Pironio and Acm 11231] 
at level 3 (that is, the maximum length of all tuples of op¬ 
erators is 3 in the sequence) to get an approximation for 



FIG. 1. The optimum measurement angles for the m = n = 3 
inequality as functions of the detector efficiency 


the matrix elements and a quite tight upper bound for the 
violation. As starting values, we have chosen three ran¬ 
dom settings on the Bloch sphere (the same for each party) 
along with a random symmetric 3-qubit state. Note the ini¬ 
tial qubit measurements do not confine in general the state 
to a 3-qubit state during the iteration procedure. We have 
repeated the above iterative procedure many times with dif¬ 
ferent (randomly generated) starting values. The best in¬ 
equality we have got is the one given above. We note that 
the only assumption we have made is that we allowed no 
single-party marginals, we actually allowed non-symmetric 
inequalities, and the procedure constrained neither the state, 
nor the measurement operators, not even the dimensionality 
of the Hilbert space. 



FIG. 2. The maximum violation of the m = n = 3 inequality as 
a function of the detector efficiency 

We applied the see-saw algorithm ll24ll to the inequal- 
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ity using three-qubit states. From this calculation it turned 
out that this space was enough to get the violation and 
'Hcrit = 1/2. We could also hnd out from this calculation 
the structure of the state and the behavior of the measure¬ 
ment operators. We have also derived a Bell inequality with 
TO = 3 and n = 4, and we have studied its behavior with the 
see-saw algorithms. These results, and the ones for to = 2 
and n > 2 by Larsson and Semitecolos in made it possi¬ 
ble for us to make the generalization and to get the family 
of inequalities presented in this paper. 

We have calculated numerically the optimum measure¬ 
ment angles and the maximum violation as functions of the 
detector efficiency above the threshold efficiency pcrit = 
0.5 for the inequality given in Eq. (l39l l. To get more ac¬ 
curate numbers, we have not used the see-saw algorithm to 
get these results, only the knowledge we had gained from it 
about the form of the solution. We have used the analytical 
expressions given by Eqs. (fTTT i and (l22l i for the matrix ele¬ 
ments, and optimized numerically the measurement angles 
and the mixing angle a. 

The results are shown in Eigs. [T]and|2] The hrst mea¬ 
surement angle tends to zero at the fastest pace, it behaves 
as the square of p — rjcrit- If we choose its value exactly 
zero, the violation changes very little. The second measure¬ 
ment angle is proportional with p — rjcrit, while the third 
one behaves as the square root of p — pcrzt- The violation is 
is proportional with the sixth power of p — pcrtt- 


VII. SUMMARY 

In the present paper we have dehned a family of binary 
outcome n party Bell inequalities with m < n measurement 
settings for each party, whose to > 2 members can be vio¬ 
lated by less detection efficiency than any other inequalities 
known so far. This gives upper bounds for the minimum 
value 77 * of the critical efficiency achievable for such in¬ 
equalities. The family is the gen eralization of the one given 
by Larsson and Semitecolos lU^ for to = 2. Unfortunately, 
for TO > 2 there is still a gap between our upper bounds and 
the best lower bounds n/((n — 1)to -f 1) given by Massar 
and Pironio d. Eor TO = n = 3 our upper bound is 1 /2, 
while the lower bound is 3/7. As we have found no better 
Bell inequality while making assumptions neither about the 
symmetry nor about the violating state, p* may well agree 
with our bound for both this case, and also for to = 3; 
n = 4. Eor n = m our bound Eq. ® gives p* < 2 /(to-|- 1), 
while from Ref. in > 1 /(to — 1 -f 1 /to) follows, a 
factor of two difference for large to. Buhrman et al. in 
gives a worse upper bound of p* < 8 /to. When n 00 , 
V* ^ 2/(to -f 2), as it follows from Eq. (| 8 ll. In this case 
the lower bound is p* > 1 /to, so there is still a factor of 
two difference for large to. The upper bound of Ref. d 
behaves very similarly to ours, it is p* > 1/to (see Eq. (4) 


in their paper). It is interesting to note that 77 * does not ap¬ 
proach zero if n alone goes to inhnity, it is proportional to 
1 /to. However, if to 00 while n remains hnite, the lower 
bound of Ref. d goes to zero. Ref. in does not give a 
usefu l up per bound for this situation, while in the case of 
Ref. lU^ the critical efficiency remains hnite, proportional 
to 1/n. Unfortunately, the present work tells nothing about 
the m > n case, so even the qualitative behavior of 77 * in 
this case remains open. 

Eor the present family of Bell inequalities the violat¬ 
ing state approaches a product state near the threshold ef- 
hcien^. The same behavior has been observed by Vallone 
et al. fl] for several other Bell inequalities, e. g. for the two- 
qubit chained inequalities with any number of measurement 
settings. They even concluded that the complementary be¬ 
havior of entanglement and nonlocality might be general. In 
the multipartite case this can not be so, as the inequalities 
given in Refs. d and d are violated by the maximally 
entangled GHZ and W states near rjcrit, respectively. In the 
bipartite case, a counterexample is shown in the Appendix. 
However, it is still possible that the conjecture by Vallone 
et al. is correct for inequalities with the smallest critical ef- 
hciency, the ones with rjcrit = V* ■ 

Our multipartite Bell inequalities involve high-order cor¬ 
relators among the particles, which are challenging to ac¬ 
cess in case of many particles. However, recently exper¬ 
imentally more friendly binary-outcome Bell inequalities 
have been constructed involving only two-body correla¬ 
tors . We pose it as an interesting problem to hnd such 
low-order Bell inequalities which are suited to Bell viola¬ 
tion using detectors of low efficiencies. 
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Appendix A: Two-party Bell inequality where the maximal 
robustness of nonlocality is attained with a maximally 
entangled state 

To prove the title of the appendix, let us consider the fol¬ 
lowing Bell inequality: 

h =P(Ai, Bi) + P(Ai, B 2 ) + P{.A2, Si) - P{A2,B2) 
- cP(Ai) - cP(Bi) < 2 - 2c, (Al) 

where c = ~ 0.9136. We choose that both Al¬ 

ice and Bob always output “0” in case of non-detection, 
hence all measurements have binary outputs (it turns out 
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that other choices will not lead to lower detection effi¬ 
ciency thresholds). Then, the measurement probabilities 
are modihed according to P{Ax,By) —>■ rf‘P{Ax,By), 
P{Ax) —?> rjP{Ax) and similarly for Bob. Introducing 
these expressions in Eq. (lAlb and dividing by rp, we ob¬ 


tain the detection-efficiency-dependent inequality: 

12(77) =P(Ai, Bi) + P(Ai, B2) + P { A 2 , Bi ) - P{A2,B2) 

(A2) 

77 r] r]‘^ 

By setting rj = c, we get the Clauser-Horne inequality 
on the left-hand-side, whereas the right-hand-side becomes 
(v/2 —1)/2, which is just the maximal quantum violation of 
the Clauser-Horne inequality attainable with a pair of max¬ 
imally entangled qubits vQm . From this it follows that for 
any 77 approaching c ~ 0.9136 from above the maximal 
violation of the Bell inequality (lAlb is given by a state con¬ 
verging to the maximally entangled two-qubit state. 
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